
INKIA 



INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE 



OC 

c 

(N 



Analytical Solution for Wave Propagation in 
Stratified Poroelastic Medium, Part II: the 3D Case 



Julien Diaz — Abdelaaziz Ezziani 



in 

(N 

< 

TO 



> 

o 

OC 
C 



TO 



N° 6596 

Juillet 2008 

.Theme NUM 




o 

z 

LU 

+ 



OS 

in 

d: 






05 

to 



CO 
CO 



INSTITUT NATIONAL 

DE RECHERCHE 

EN INFORMATIQUE 

ET EN AUTOMATIQUE 




INKIA 



centre de recherche 
BORnEAUX - sun QUEST 



Analytical Solution for Wave Propagation in Stratified 
Poroelastic Medium. Part II: the 3D Case 



Julien Dia2lj[!J, Abdelaaziz Ezziani^^ * 

Theme NUM — Systemes numeriques 
Equipe-Projet Magique-3D 

Rapport de recherche n° 6596 — Juihet 2008 — [371 pages 



Abstract: We are interested in the modehng of wave propagation in poroelastic media. We 
consider the biphasic Biot's model in an infinite bilayered medium, with a plane interface. 
We adopt the Cagniard-De Hoop's technique. This report is devoted to the calculation of 
analytical solution in three dimension. 
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Solution analytique pour la propagation d'ondes en milieu 
poroelastique stratifie. Partie II : en dimension 3 

Resume : Nous nous interessons a la modelisation de la propagation d'ondes dans les milieux 
infinis bicouches poroelastiques. Nous considerons ici le modele bi-phasique de Biot. Cette 
seconde partie est consacree au calcul de la solution analytique en dimension trois a I'aide de 
la technique de Cagniard-De Hoop. 

Mots-cles : Modele de Biot, ondes poroelastiques, solution analytique, technique de Ca- 
gniard de Hoop. 
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Introduction 

Many seismic materials cannot only be considered as solid materials. They are often porous 
media, i.e. media made of a solid fully saturated with a fluid: there are solid media perforated 
by a multitude of small holes (called pores) filled with a fluid. It is in particular often the 
case of the oil reservoirs. It is clear that the analysis of results by seismic methods of the 
exploration of such media must take to account the fact that a wave being propagated in such 
a medium meets a succession of phases solid and fluid: we speak about poroelastic media, 
and the more commonly used model is the Blot's model [H El E]. 

When the wavelength is large in comparison with the size of the pores, rather than regarding 
such a medium as an heterogeneous medium, it is legitimate to use, at least locally, the theory 
of homogenization [11[13]. This leads to the Blot's model [HEllS] which involves as unknown 
not only the displacement field in the solid but also the displacement field in the fluid. The 
principal characteristic of this model is that in addition to the classical P and S waves in a 
solid one observes a P "slow" wave, which we could also call a "fluid" wave: the denomination 
"slow wave" refers to the fact that in practical applications, it is slower (and probably much 
slower) than the other two waves. 

The computation of analytical solutions for wave propagation in poroelastic media is of high 
importance for the validation of numerical computational codes or for a better understanding 
of the reflexion/transmission properties of the media. Cagniard-de Hoop method [Sj [7] is a 
useful tool to obtain such solutions and permits to compute each type of waves (P wave, S 
wave, head wave...) independently. Although it was originally dedicated to the solution to 
elastodynamic wave propagation, it can be applied to any transient wave propagation prob- 
lem in stratified medium. However, as far as we know, few works have been dedicated to the 
application of this method to poroelastic medium, especially in three dimensions. 

In order to validate computational codes of wave propagation in poroelastic media, we 
have implemented the codes Gar6more 2D [11] and Gar6more 3D ^2j which provide the com- 
plete solution (reflected and transmitted waves) of the propagation of wave in stratified 2D 
or 3D media composed of acoustic/acoustic, acoustic/elastic, acoustic/poroelastic or poroe- 
lastic/poroelastic layers. The codes are freely downloadable at 



http : //www . spice-rtn . org/library/sof tware/Garciinore2D 
and 



http : //www . spice-rtn . org/library/sof tware/Gar6more3D[ 



We will focus in this paper on the 3D poroelastic case, the two dimensional and the acous- 
tic/poroelastic cases are detailed in [8l|9l[10]. The outline of the paper is as follows: we first 
present the model problem we want to solve and derive the Green problem from it (section 
1). Then we present the analytical solution to the wave propagation problem in a stratified 
2D medium composed of an acoustic and a poroelastic layer (section 2). Finally we illustrate 
our results through numerical applications (section 3). 
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1 The model problem 

We consider an infinite two dimensional medium (Q = R ) composed of two liomogeneous 
poroelastic layers Q^ = R x] — c>o,0] and Q~ = R x [0,+oo[ separated by an horizontal 
interface T (see Fig. [1]). We first describe the equations in the two layers ( §l.ip and the 
transmission conditions on the interface F ( §1.2|) . then we present the Green problem from 
which we compute the analytical solution ( §1.30 . 



First Layer 

0+ 



z = 



Second Layer 



Figure 1: Configuration of the study 



1.1 Poroelastic equations 

We consider the second-order formulation of the poroelastic equations [H EJ [3] : 

in nx]0,T], 

in nx]0,T], 

in nx]0,T], 

P + f3V -Us + V-W = F„, mnx]0,T], 

m 

Usix, 0) = 0, W{x, 0) = 0, in Q, 

tlsix, 0) = 0, W{x, 0) = 0, in Q, 



Pftjs + p^nW + ^W + VP = F^., 

S = AV ■ [/, Is + 2pe{Us) -f3Pl3, 
1 



with 



(1) 



dT,- ■ 
(V • S)j = y^ ''' y i = 1,3. As usual /s is the identity matrix of ^^2(2^)) 



and s{U s) is the solid strain tensor defined by: 
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In ([T]), the unknowns are: 

• Us the displacement field of solid particles; 



W = 4>{U f — Us), the relative displacement, Uj being the displacement field of fluid 
particles and (p the porosity; 



• P, the fluid pressure; 

• S, the solid stress tensor. 

The parameters describing the physical properties of the medium are given by: 

• p = (j) pf + {1 — (t>)ps is the overall density of the saturated medium, with ps the density 
of the solid and pj the density of the fluid; 

• Pw = o-pf/cj), where a is the tortuosity of the solid matrix; 

• /C = k/i], where k. is the permeability of the solid matrix and r/ is the viscosity of the 
fluid; 



m and P are positive physical coefficients: P = 1 — K^/Ks 

and m = [(p/Kf + (/3 — (l))/Ks]~ , where Ks is the bulk modulus of the solid, Kf is the 

bulk modulus of the fluid and Ki, is the frame bulk modulus; 



• /x is the frame shear modulus, and X = Kb — 2p/3 is the Lame constant. 

• Fu, Fw and Fp are the force densities. 

To simplify this study, we consider only the case of a compression source 

F„(x, y, t) = fuV{5cc Sy 6z-h) f{t) and F^(x, y, t) = U^{6^ Sy d^^h) fit) 

and a pressure source Fp = fp6x6y6z-h f{t), where /«, fw and fp are constant and / is a 
regular source function in time. We can generalize this approach for other types of punctual 
sources such as for instance 

Fu = fvS X (4 5y 5^-h v) f{t) and Fy, = f^jV x {6^ 6y 5^-h v) f{t) 

where t? is a vector in IR . 
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1.2 Transmission conditions 

Let n be the unitary normal vector of T outwardly directed to 0,^ . The transmission condi- 
tions on the interface T between the two poroelastic medium are [6j: 

r ut = u-, 

W+ ■ n = W ■ n, 
P+ = P-, 



(2) 



1.3 The Green problem 

We won't compute directly the solution to ([2]) but the solution to the following Green problem: 



1 






m 



± 



p^ + P^V-uf + V-w^ = fp 5x 6y 5,.h fit), 



in n^x]0,T], 
in O=^x]0,r], 
mn^x]0,T], 



m7 = u+, onrx]0,r] 

w n = w^ n, onTxJOjT] 

J9"=p+, onrx]0,T] 

a-n = a^n, onrx]0,r]. 



(3a) 

(3b) 
(3c) 
(3d) 

(3e) 
(3f) 

(3g) 
(3h) 



The solution to ([T]) is then computed from the solution to the Green Problem thanks to a 
convolution by the source function. For instance we have: 

P^ix,y,t) =p^{x,y,.)* f{.) = / p+{x,y,T)f{t-T)dT 

Jo 

(we have similar relations for the other unknowns). We also suppose that the poroelastic 
medium is non dissipative, i.e the viscosity rj^ = 0. Using the equations (|3c|3d|) we can 
eliminate a and p in ([3]) and we obtain the equivalent system: 

' p^uf + pf w^ - a^ V(V • u^) + ^=^ V X (V X uf) - m^l3^V{V ■ w^) 

= [fu - P^m+fp)V{6x 6y 6,-h) St, (4) 

t pfut + p± ib^ - m^/3^ V(V • uf) - m± V(V • w^) = {U - m+fp)V{6x 6y 5,_0 6t, 
with a~ = A^ + 2fi^ + m~/3~ . 
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And the transmission conditions on T are rewritten as: 



^sx ^sxy 



"'sy "'syi 



W, = W' 






(5a) 
(5b) 
(5c) 
(5d) 
(5e) 

(5f) 

(5g) 

(5h) 



l^^idzuty + dyut^) = fi {dzu,y + dyu,^), 

(A- + m+/3+^)V • w+ + 2/2+9^^+ + m+/3+V • w+ = 

(A~ + m~ (3~ )V • w7 + 2/i~92u7z + m^ (3^V ■ w . 

We split the displacement fields uf and w^ into irrotational and isovolumic fields (P-wave 
and S-wave): 



u 



± 



ve^ + V X ^± ; w^ = ve^ + V X * 



(6) 



The vectors ^„ and ^^ are not uniquely defined since: 

Vx (*^ + VC) = Vx *^, V£g{u,u;} 
for all scalar field C . To define a unique ^^ we impose the gauge condition: 

V • *^ = 0. 
The vectorial space of ^ ^ verifying this last condition is written as: 



* 



dy 




r 92, 1 


-dx 


^?i + 


dl 







-92 -92 

L '-'xx ^yy -1 



^ 



£,2' 



where ^^ ^ and ^^ 2 ^^^ two scalar fields. The displacement fields uf and tu^ are written in 
the form: 



92 



ut 


= v@t + 


92 

^yz 

-92 -92 
'^xx ^yy 

r di 


w"^ 


= ve± + 


dlz 

-92 -92 






* 



UI,1 



9j, 


dy 

-dx 




A^^,2 



(7) 



A^^,2- 
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We can then rewrite system (j3]) in the following form: 

A-Q- - B-Ae- = 0, in f^" x]0, T] 

K^-^t^K^ = ^^ ie{i,2}, inO±x]o,r] 

in n^x]0,T] 



(8) 






Pw 



where G± = (e-,G-)*,F 
matrices: 



ifu - /?+m+/„, /, 



p^ JW 



m 



~fpY, A and i? are 2x2 symmetric 



A* 



Pf Pt 



B' 









m^ 



and 



vs 



1 



Atp^ 



± ± ±2 



is the S-wave velocity. 



We multiply the first (resp. the second) equation of system ([8]) by the inverse of A~^ (resp. 
A-). The matrix A+~^B+ (resp. A-'^B') is diagonalizable: A^~^B^ = V^D^V^~^, 
where V^ is the change-of-coordinates matrix, D^ = diagiVpr , Vp^ ) is the diagonal matrix 
similar to A^ B^, Vpr and Vp^ are respectively the fast P-wave velocity and the slow P- 
wave velocity (Vp^ < Vpr). 



^±-i«± 



Using the change of variables: 

we obtain the uncoupled system on fast P-waves, slow P-waves and S-waves: 

' ^+ - D+A^+ = 4 6y 6^-h St F+, in rj+x]0, T] 
^- - D-A<p- = 0, inO-x]0,T] 

Ki-^f^Ki = ^^ ^G{1,2}, in^±x]0,T] 



(9) 



(10) 



Pw 



in n^x]0,T] 



with F+ = {A+V+)-^F = (F+^,F+J*. 

Using the transmission conditions (f5al) - ()5bp . (|5fj)-([5gl) and the change of variables ([7j), we 
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obtain: 



d,e+ + aiV'+i - dyiA^:,) = d,e- + a^c, - dyiA^p-,), on r 



flla) 



dyet + d'y,^+ + 5,.(AV+ ) = dyQ- + ^'y,^p-, + a,(AV-o), on r, (lib) 



Ai" 



29^9+ + 5.(aL - A^)^l, - d^A^X,) 



fi- 29^9- + d.idl - A^)^-, - d^A^Za) ' 



on r, (lie) 



^.+ 29^,9+ + a,(aL - Ax)^+i + aiAv&+ 



onr, 



fi- 29,^,9- + 5,(ai, - Ax)*-, + 94 AvI/-. 



(lid) 



with A^ = d'^^ + (9^ . Applying the derivative dy to the equation pia|) (resp. picp l. dx to 
the equation (lllbp (resp. (jlldp ) and subtracting the first (resp. the third) obtained equation 
from the second (resp. the fourth) one, we get: 

Ax ( AV'+2) = A± (AV-2) , on T, (12a) 

^+(a,Ax)A*+2 = ^-(a,Ax)AM/-2, on T, (12b) 

moreover, using the third equation of (llOp . we have *„ 2 satisfies the wave equation: 

^n,2 - ^s^'A*i2 = 0, inO±x]0,T] 
and, since u^ and w^ satisfy, at t = 0, uf- = iig = w^ = w = 0, we obtain: 

^±2 = 0, inl7±x]0,T], (13) 

and from piap -( ITTbll we deduce the transmission condition equivalent to (I5ap and ()5bp : 

d.&t + dliP+=d.e- + dli:-„ onr. (14) 



In the same way, using the equality ()13p . we can show that the two transmission conditions 
(jllcp and (jlldp are equivalent, which gives us: 

M+ (2^0;^ + dxidl - Ax)*+i) = A^- (2dle- + a,(aL - Ax)*-i) , on r. (15) 

We can then reduce the transmission conditions ([5]) to 6 equations: (fTH [5ct I5d|5e|15|5h|) . 

Finally, we obtain the Green problem equivalent to ([3]): 

^t-V+^A<^+ = 6x6y6,-h6tF+, t^{Pf,Ps} z>0 
^+ - F<t^A$+ = z > 

^r - Vr'^A^r = o, i G {P/, Ps, S} z <0 

[ mtp ^Ps^ ^S^ ^~Pf^ ^~Ps^ ^~S)=0^ ^ = 
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where we have set ^g = ^^ ^ in order to have similar notations for the Pf, Ps and S waves. 
The operator B represents the transmission conditions on T: 



B 



(Hf\ 




<^t>s 




^1 




^p/ 




^Ps 




V ^5 / 





vt, d. 


Vt2d. 


dlz 


-Pn9. 


-^12 ^^ 


-C 


K, d. 


Vt2dz 


-A^ 


-^fi^. 


-^1-2 5. 


A^ 


^41 


Vt^dz 
^42 


Pw 




-^2"i 5^ 

^44 


— / 22 ^z 
S45 





^51 


^52 


^53 


S54 


^55 


^56 


^61 


^62 


-2/i+9,(Ax) 


^64 


^65 


2/x-9,(Aj_ 





\^h] 




^ts 




^i 




^-pf 




^Ps 




.^5 . 



where P--, i,j = 1,2 are the components of the change-of-coordinates matrix V^ and 



B, 



41 



m+(/3+7^++P+).o2 . „._ _ m+(/3+7'+ +7^2+2) ^2 



V^i 



■5/^; ^42 



P/ 



V, 



d, 



tt^ 



Ps 



B44 



"^"(r^ri+^2"i)«2 . 



ai ; ^4 



"i-(r^r2+^2"2);.2. 



^p7 ^i . 

^51 = 2A.+P+ 9i ; B,2 = 2/i+P+ 9^ ; B53 = /U+5,(C " A±); 
S54 = -2//-pn ai ; B55 = -2fi-V^2 ^L ; ^56 = -^"5,(9^ - A^); 
(A+ + m+P+^)Vti + m+ZJ+P^l .2 



^61 
^62 
i364 
^65 






-5/, + 2/.+p+a 



2 . 

zz^ 



(A+ + m+/3+ )7'i+2 + m+f3+V. 



22 q2 






5,l + 2^.+p+25..; 






F, 



Pf 



(A- + m-/3- )Pr2 ±rn2£r^ ^2 + 2/.-pr2C. 



T^, 



Ps 



To obtain this operator we have used the reduced transmission conditions (I141 [5c ll5d|5e|15|5hp . 
the change of variables (|6I9|) and the uncoupled system (|1U|) . 

Moreover, from the unknowns ^pf, ^p^ and $^ we can determine the solid displacement 
uf- and the relative displacement w^ by using the change of variables presented below. 



2 Expression of the analytical solution 

Since the problem is invariant by a rotation around the z-axis, we will only consider the case 
y = and x > 0, so that the y-component of all the displacements are zero. The solution for 
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11 



y 7^ or X < is deduced from the solution for y = by the relations 

X 



Usxix, y,z,t) 



-uf^{y/x^~+y^,0,z,t) 



uty{x,y,z,t) = uf^{^Jx'^+y'^,{),z,t) 

y x^ + y^ 



U7.{x,y,z,t) = uf^ ( \/x2 + y^, 0,z,t) 



^sz\ 



(17) 
(18) 
(19) 



To state our results, we need the following notations and definitions: 

1. Definition of the complex square root. For q £ C\IR~, we use the following 
definition of the square root g[q) = q^'"^: 

g{qf = q and ^e[g{q)] > 0. 

The branch cut of g{q) in the complex plane will thus be the half-line defined by {q £ 
]R~} (see Fig. [2]). In the following, we use the abuse of notation g(q) = i\/—q for 
qeJR-. 



Qm{q) > 


^ 


+7T^ 


^ 




' ^ 


"V 


J 3f?e( 



Figure 2: Definition of the function x i-^ (x)^' ^ 



2. Definition of the fictitious velocities For a given <; G M, we define the fictitious 
velocities Vj (g) for i G {P/, Ps^ S\ by 



Vf := Vf (g) = Vt 



1 



1 + K- 



±2^2' 



These fictitious velocities will be helpful to turn the 3D-problem into the sum of 2D- 
problems indexed by the variable q. Note that V^ (0) correspond to the real velocities 



± 



3. Definition of the functions k^ . For i G {P/, Ps, 5*} and (g^, g^) G C x IR, we define 
the functions 



H :=^ {(lx,%) 



V: 



^2+11 + 11 



1/2 



1 



1/2 



yt\%) 



+<il 
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4. Definition of the refiection and transmission coefficients. For a given q = 
{qx,qy) G C X IR, we denote by TZpfpf{q), TZpfPsiq), Tlpfsiq), '^PfPfiq), '^PfPsiq), 
and Tpfs{q) the solution to the linear system 



A{q) 



T^pfpfiq) 




T^pfPsiq) 




T^pfsiq) 


1 


'Tpfpfiq) 


24/g)F+/ 


'^pfPsi.q) 




'Tpfsiq) 





-4/q)P- 



11 



m 



Vpf 



-Kpf{q)V^i 



2iga;/i+K+^(q)P 



11 



(A+ + m+/3+^)P+ + m+(3+V^^ 



V, 



2^+P+4/(g) 



Pf 



and by TZpsp/iq), TZpsPsiq), TZpssiq), TpsPf{q), TpsPs{q) and Tpss{q) the solution to 
the linear system 



A{q) 



T^PsPfiq) 

TlpsPsiq) 

T^Pssiq) 

'TpsPfiq) 

TpsPs{q) 

Tpssiq) 



24,{q)V^: 



m ' 



V, 



^qxV^2 

-^Psil)'Pt2 
-4si<l)T'22 

m+rt^ + r, 



'22) 



Ps 



2ig^/i+K^^(q)P- 



12 



(A+ + m+/3+^)P+ + m+[3+V. 



22 



V, 



2fi+V+K+/{q) 



Ps 



where the matrix A{q) is defined for q = (qx, qy) G C x IR by: 



Mq) 



-'^QxVii 



-«P/(9)^^i 
^41 (q) 

Ai(q) 
Ai(q) 



-i qxV{2 



-4si<i)'Pi2 

-^Psil)'Pt2 
^42 (q) 

A2(q) 

A2(g) 



^QxUgiq) 

Ql + Qy 



" Pw 





A3(q) 
Am{q) 



i 9x^11 



-K'Pf{q)'Pii 

-Kpj{q)V2i 

^44(9) 
A4(g) 
A4(q) 



i QxV^2 



-^Ps(.<l)'P: 



12 



-'^;.(q)^2"2 iQl + QlV:^ 



^QxUgiq) 



-<& - qI 



2-.Pj_ 



AA5{q) 
A55{q) 
Ae^iq) 





A56{q) 
Amiq) 
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with 

TTl Tfi 

-^4i(g) = —ill^^'Pti+'Pti]-^ -^44(9) = -— T[/3"^ri + ^2"i]; 
pf pf 

TTl TTl 

^42 (q) = ^72 [/3'^^l'2+^2"'2]; -^45(9) = —[l^"Pl2 + P22\-^ 

Ai(q) = 2ig^^+K+^(q)P^+i; ^54(9) = 2ig^^~Kp^(g)P{"^; 

^2(9) = 219^^+^+^(9)7^+; A5(q) = 2\q^^x- K^^{q)P{^] 

A3(q) = -iM"^9x(K5 (q) +g^ + 9^); A6(q) = i^"9x(K5 (q) + g^ + g^); 

, r„^ - (A+ + m+/3+^)P+ + m+/3+P+ 2 + . 

>l62(qj - :p2 ^ ^/^ '*Psl9J n2' 

^Ps 

A3(Q) = -2((?^ + (?^)^+4(q); 

A4(q) = 7^^ — - 2fi Kpj (q)Pii; 

^Pf 

M^{q) = -^ - 2/U Kp^ (g)Pi2; 

^Ps 

We also denote by T^nax the greatest velocity In the medium: 

V^max = max(F/^, Vj^^, V^, Vpj, Vp^, Vg). 
We can now present the expression of the solution to the Green Problem: 
Theorem 2.1. The solid displacement in the top medium is given by 

uf{x,0,Z,t) = Upj{x,Z,t) + Upj:pj:{x,Z,t) + Upjp^{x,Z,t) + Upjg{x,Z,t) 

+ u^^{x,z,t) + u^^p^{x,z,t) + u^gpg{x,z,t) + u^_,g{x,z,t) 
and the solid displacement in the bottom, medium is given by 

U~{x,0,Z,t) = U^j:pj:{x,Z,t) + U'^jp^{x,Z,t) + U^j:g{x,Z,t) 

+ U+^pj:{x,Z,t) + u'^^p^{x,z,t) + uj,^g{x,z,t), 
where 
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Upr is the solid displacement of the incident Pf wave and satisfies: 

rtiF+f xtH{t - to) 



u 



Pf. 



^{x,z,t) :-- 



V, 



2 47rr3 



uj,j:Jx,z,t) :-- 



Pf 

'PiiFpf {z - h)tH{t - to) 
T/+ ^ 47rr^ 



where H denotes the usual Heaviside function. Moreover we set r = (x^ + (z — h)"^)^'^ 
and to = r/Vpr denotes the time arrival of the incident Pf wave at point {x,0,z). 



• Mp^, is the solid displacement of the incident Ps wave and satisfies: 

P+^+, xtH{t - to) 



UpsJx,z,t) :-- 



)+ Z7+ 



u^^ {x,z,t) :-- 



Vt^F^, {z - h)tH{t - to) 



V, 



Ps 



47rr^ 



We set here r = (x^ + (z — /i)^)-^'^ and to = r/Vp^ denotes the time arrival of the 
incident Ps wave at point (x, 0, z). 



u'^rpr is the solid displacement of the reflected PfPf wave (the Pf reflected wave 
generated by the Pf incident wave) and satisfies: 



''PfPf,: 



^{x,z,t) = 'PtiFtf j^ 

r 

'^PfPf,z\^-''^i'^) ~ F'uFpf 

ift^^ <t<to and-> -f^, 

' •'max 

'^PfPf,x\''^i^i) ^ 



<ji(t) 9m 



v{t, q)KJ,r{v{t, q))npfpf{v{t, q)) 



qi{t) Qm 



7r'^ry/q'^+q^{t) 
i^Pf {'^{t,Q))'^PfPf{'^{t,Q)) 



dq, 



7r2rvV + g!(i) 



■dq, 






gi(t) 9m 



v{t, q)K+r{v{t, q))TZpfPf{v{t, q)) 



qo(t) 



vr- 



'ry/q"^ - q^{t) 



■dq 



'^PfPf.zy^^i ^1^) 



ftlPps f 



(t) ^e 



7(t, q)Kpf{j{t, q))'JZpfPf{7{t, q)) 



TT^r^MW 



■dq, 



qi(t) Qm 



i^Pf {v(t,q))npfpf{v{t,q)) 



qoit) Tr^r^/q'^-q^{t) 



(t) SRe 



i^Pf {i{t,q))T^pfPf{i{t,q)) 



TT^ry^q^{t) - q^ 



dq 
dq, 
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V, 



ifto<t< tu and - > 77 



Pf 






90 



(t) SRe 



lit, (l)i^'l>f{l{t, q))T^PfPf{l(t, q)) 



u 



+ 

PfPf,z 



{x,z,t) = -Vt^F+^ 



<?0 



(t) 'Re 



7r'^r^/q^{t) - q^ 
i^Pf i7it,q))'R-PfPfi7it,q)) 



■dq, 



IT- 



■Wqm-Q' 



dq, 



X ''^Pf X ^Pf 

if tj^ < t and — > — or if to < t and — < — and 

f •'max ^ ''max 

upfpf{x, z,t) = else . 

We set here r = (x^ + (z + h)'^)^''^ and to = r/Vpr denotes the arrival time of the 
reflected PfPf volume wave at point (x,0, z), 



1 1 \x 

th, = {z + h)J^^-^^:^ + 



T/+ 2 y2 Y 

V p f max * max 



(20) 



denotes the arrival time of the reflected PfPf head-wave at point {x, 0, z) and 

tho 



1 



1 



Z + h^V+/ ^rLx 



(21) 



denotes the time after which there is no longer head wave at point (x,0, z), (contrary to 
the 2D case, this time does not coincide with the arrival time of the volume wave). We 
also define the functions 7, v, qo and qi by 



xt z + h If^ 1 

7 : {t G iR 1 1 > to} X iR H^ C := 7(t, gj;) = i — + 



^ V^' ^t){%) 



Pf\Hy; 



V : {t e IR\th-^ <t < th2}xlR^ C := v{t,qy) = -i 



z + h 



Vffiqy) ^ 



t^ X 
+ ^t 



qo : m ^ m := qo{t) 



\ 



t2 1 



r^ Y 



Pf 



and 



qi : m ^ m := qi{t) 



\ 



t-iz + h)^ 



1 1 



T/+2 y2 I y2 • 

Vpf "^ max ' •'tv,,^ 



u^fp^ is the solid displacement of the reflected PfPs wave and satisfies: 



RR n° 659 



u 



PfPs,x 



[x,z,t) 



Vt2Fpf /■'^iW 



vr^ 



Re 



dv 

i v{t, q)TZpfPs {v{t, q))-Q^{t, q) 



dq, 



'^PfPs,z\^^^^^) 



vt^F^f r^iw 



vr^ 



Re 



i^Psi^it, q))'^pfPs{v{t, q))-Q^(.t, q) 



dq, 
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ifth^ <t<to and \Qm[y{to,0)]\ < 



V ' 

* max 



'^PfPs,x\^^^^^) 



vt2Fpf r*(*) 



TT^ 



3f?e 



vt^F^f r'/iW 



TT^ 



3f?e 



90 (t) 



'^PfPs.z\^-'^^^) 



Vt2Fpf /■*(*) 



TT^ 



3f?e 



7?+^+, f'/iCt) 



12^ p/ 

7r2 



ifto<t<th2 and \Q'm[y{tQ,0)]\ < 



^e 



90 (t) 
1 



(?7 
i-i{t,q)TZpfPs{l{t,q)) — {t,q) dq 

dv 
iv{t,q)TlpfPs{v{t,q)) — {t,q) dq, 






Kpsi^i'^^ q))T^pfPs(.v{t, q))-Q^(.t, q) 



dq 
dq, 



V 



U 



PfPs,x 



U 



+ 
PfPs,z 



[x,z,t) 



[x,z,t) 



p+F+^ r9o(t) 



vr^ 



3f?e 



T'uFpf /"'"(*) 



vr^ 



3f?e 



O7 
i-f{t,q)TZpfPs{-i{t,q)) — {t,q) dq, 



dq. 



ifth2 < t and |9m [7(^0, 0)]| < — — or if to < t and \Qm[j{to,0)]\ > — — 

''max ^max 

and Uprp^{x, z,t) = else. 

to denotes here the arrival time of the reflected PfPs volume wave at point {x, 0, z) 
(its calculation is similar to the calculation of the arrival time of the transmitted wave, 
see the appendix of I9j), 



th^ = h 



1 



y+/ K^ax 



+ z. — 



1 \x\ 

+ 



Vp^ ''max ''max 



(22) 



denotes the arrival time of the reflected PfPs head wave at point {x,0,z), 



th2 



K' + z' + hz[^f^+f^]+x' 

_h 1 z_ 

ci C2 



(23) 



denotes the time after which there is no longer head wave at point {x,0,z), where 



1 1 



ci 



V^/ ^^ax 



and C2 



1 1 



Vn„ "mux 



' Ps 



The function qo : [to J +00] ^^ IR^ is the reciprocal function of to '■ IR^ h^: [to, +00], 
where to{q) is the arrival time at point [x, 0, z) of the fictitious transmitted PfPs volume 
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wave, propagating at a velocity Vpr{q) from the source to the interface and at velocity 
Vpg{q) from the interface to point {x,0,z) (we refer again to I9j for details on its cal- 
culation). 

The function qi : [ti ; to] ^^ ^^ is defined by 



91 (t) 



\ 



1 1 

^ I t- z 






' -k.U 



V r, c max 



V^. 



The function 7 : {(t, (?) G IR^ x iR^ 1 1 > ^0(9)} ^^ C is implicitly defined as the only 
root of the function 

whose real part is positive. 

The function v : Ei D E2 ^^ C is implicitly defined as the only root of the function 



^'"■'■" = H^ 



1/2 



+ v^ 



h 



1 



V^/(cz) 



+ V ] + ivx — t 



such that 9m [dtv{t, q)] < 0, with 

El = {{t,q) &IR^ xIR^\th^<t<tQ andO < q < qo{t)} 

and 

E2 = {{t-,q) G ffi+ X iR+ |to < t < thi and qo{t) < q < gi(t)} 

itpr^ is the solid displacement of the reflected PfS wave and satisfies: 



^P/S,x(^'^'*) 



F+f Mt) 



7r2 / 
^ JO 



3f?e 



5f 



iv{t,q)K^{v{t,q))npfs{v{t,q)) — {t,q) 



dt 



dq, 



'"'PfS,z\^^^'^) ~ ^2 



F+ /-QiW 



TT^ 



3f?e 



dv 



{v'{t, q) + g")7^p/s(^^(i, 9))^(*' 9) 



dg, 



«/i/ii < t < to and 1 9m [7(^0,0)] I < 



^P/S,x(^'^'*) 



+ 



RR n° 



^PfS.zy^^ ^' ) 



5596 



+ 



U!";j 


" V ' 

^max 


Pf 

7r2 


Vo 


T+ 
Pf 

7r2 


Jqo{t) 


Pf 

7r2 


Vo 




/ 3f?e 



l7(t,(?)4(7(t,g))7^p/5(7(i,9))|^(t,9) 



yv{t, q)K'^ {v{t, q))npfs{v{t, q)) — {t, q) 



dq 

dq, 



TT^ 



'Qo(t) 



(z;2(i, g) + q^)TZpfsivit, <?))-^(t, ?) 



(ig 
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iftQ<t<th2 and \Qm[y{tQ,0)]\ < 



K. 



u 



PfS,x 



'x,z,t) 



-7r2 , 



qo{t) 



^e 



n{t, g)4(7(t, q))npfsh{t, q))^{t, q) 



dt 



dq, 



u 



+ 

PfS,z 



{x,z,t) 



F^f fioit) 



-7r2 , 



^e 



dj 



iY{t,q) + q')npfsi7{t,q))-^it,q) 



dq, 



if th2 < t and \Qm [7(^0, 0)]| < — or ifto<t and {Qtu [7(^0) 0)]| > 



K rn 



Kr 



an 



d Uprg{x, z,t) = else. 



to denotes here the arrival time of the reflected PfS volume wave at point {x, 0, z) 
(its calculation is similar to the calculation of the arrival time of the transmitted wave, 
see the appendix of fB/), 



thi 



h. 



V+; ^r^ax 



+ zj — 



1 \x\ 

+ 



Vo "^max ''max 



(24) 



denotes the arrival time of the reflected PfS head wave at point (x, 0, z) 



th2 



h^ + z'^ + hzi^ + ^]+x^ 



h 1 z_ 

ci C2 



(25) 



denotes the time after which there is no longer head wave at point (a;, 0,2), where 



ci 



v^; vL^ 



and C2 



T/+2 y2 ■ 

V Q ' max 



The function qq : [^0 ! +00] 1-^ SV~ is the reciprocal function of to : IR'^ 1-^: [to, +00], 
where to{q) is the arrival time at point (x, 0, z) of the fictitious transmitted PfS volume 
wave, propagating at a velocity VpAq) from the source to the interface and at velocity 
Vg{q) from the interface to point {x,0,z) (we refer again to fS^ for details on its calcu- 
lation). 

The function qi : [ti ; to] ^^ ^^ is defined by 



qi{t) 



\ 



^^1* "yvf K^a 



V p r ' ma 



y2 



The function 7 : {(t, g) G IR^ x IR^ \ t > io{q)} 1— > C is implicitly defined as the only 
root of the function 



^ii,q,t) = z 



yt\^) 



1/2 



+7^ 



+/i 



1 



1/2 



+ 7^ 



+ i'yx — t 
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whose real part is positive. 

The function v : Ei U E2 ^^ C is implicitly defined as the only root of the function 



J^{v,q,t) = z 



1 



yt\^) 



\ 1/2 

+ v'^\ +h 



1 



V^fil) 



+ V \ + ivx — t 



such that Qm [(9if (t, g)] < 0, with 

El = {{t,q) € 1R+ X lR+\th^<t<to andO < q < qo{t)} 

and 

E2 = {it,q) £ R^ X M^ \to < t < th^ and qo{t) < q < qiit)} 

'^'psPf ^"^ ^^^ solid displacement of the reflected PsPf wave and satisfies: 



''PsPf. 



xix^Z,t) 



-ptiF^s r^'^ 



vr^ 



^e 



dv 

i v{t, q)TZpsPf{v{t, Q))-g^{t, q) 



dq, 



u+^p.ix,z,t) = ii^ / ^e 



7T^ 



01) 

i^Pfi.'^ii, q))T^PsPf{v{t, q))-^{t 



dt 



dq, 



ifth^ <t <to and 1 9m [7(^0,0)] | < 



u 



PsPf,x 



'x,z,t) 



u 



PsPf,z 



'x,z,t) 



if to < t < th2 and 1 9m [7(^0,0)]! < 



,o)]|< 


1 

►^max 


KiF^s 


/■9o(t) 


7r2 


I 


niF^s 


rii{t) 


7r2 


Jqoit) 


KiF^s 


fgo{t) 


7r2 


L 


KFk 


rii{t) 


7r2 


Jqoit) 


nMI ^ 


1 



^e 



^e 



^e 



^e 



i-/{t,q)npsPf{-f{t,q)) — {t,q) dq 

dv 
iv{t,q)TZpsPf{v{t,q)) — {t,q) dq, 

I ^7 

i^Pfiiit, q))T^PsPf{i{t, 9))-^(i, q) 
KJ,f{v(t, q))npsPf{v{t, q))-Q^{t, q) 



dq 
dq, 



V ' 



U 



PsPf,x 



-p+7?+ rioit) 



u 



PsPf,z 



T^ Jo 

'a:,z,t) = -^^4^ I ^e 



^7 
il{t,q)TZpsPf{-i{t,q)) — {t,q) dq, 

I ^7 

i^Pfhit, q))T^PsPf{i{t, 9))^(i, q) 



dq. 



if th2 < t and |9m [7(^0, 0)]| < — — or if to < t and |9m [7(^0; 0)]| > 



Vrri 



v„ 



and Up^^pAx, z,t) = else. 
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to denotes here the arrival time of the reflected PsPf volume wave at point (x, 0, z), 



y V p max W Vp r ' max ' max 

denotes the arrival time of the reflected PsPf head wave at point (x,0,z), 

*^. = n_\\ ^ (27) 

Cl ~'~ C2 

denotes the time after which there is no longer head wave at point {x,0,z), where 



Cl = , / ^ — 7-5 — and C2 



,y+/ K^ax ^ yy+Z K^ax' 

T/ie function go ^ [*o ; +00] ^^ -K"*" is the reciprocal function of to : IR^ 1-^: [to, +00], 
where to{q) is the arrival time at point {x, 0, z) of the fictitious transmitted PsPf volume 
wave, propagating at a velocity Vpg{q) from the source to the interface and at velocity 
Vpf{q) from the interface to point {x,0,z). 

The function qi : [ti ; to] ^^ ^^ ^-5 defined by 



Qiit) 



-^ It- z 



\x^y yy+/ i^^ax v^p/ ^--/ ^--' 



The function 7 : {{t,q) £ M^ x IR^ \ t > to{q)} h^ C is implicitly defined as the only 
root of the function 



1/2 / \ 1/2 



1 o 1 . I i 2 



J'{-f,q,t) = z\ +Y +h\ ^2. +7 +i7X-t 



whose real part is positive. 

The function v : EiU E2 ^^ C is implicitly defined as the only root of the function 



1/2 / \ 1/2 



1 o\ . / 1 2 



J^(v,q,t) = z\ K \-v +h\ ^ \- V +ivx — t 



such that Qm[dtv{t,q)] < 0, with 

El = {{t,q) G iR+ X iR+|tftj < t < to andO < q < qo{t)} 

and 

E2 = {(t, g) G iR+ X iR+ I to < t < t/j^ and qo{t) < q < gi(t)} • 
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ut> p is the solid displacement of the reflected PsPs wave and satisfies: 



Jo 

Jo 
ift^^ <t<to and->-^, 

' •'max 



qi(t) 9m 



iv(t, q)K'^^{v{t, q))npsPs{v{t, q)) 



qi(t) Qm 



7T^r^/q^ + q^{t) 
>ips {v{t,q))TZpsPs{v{t,q)) 



dq, 



dq, 



^PsPs.yV^'^'*) 



' 12^ Ps 



' 12-^ Ps 



' 12^ Ps 



' 12^ Ps 



qi{t) Qm 



TrVvV + ^IW 



[iv{t, q)KJ,^{v{t, q))TZpsPs{v{t, q)) 



90 



(t) Jfte 



T^2^^q2 -(git) 

i-/{t, q)KJ,^{j{t, q))npsPs{l{t, q)) 



dq 





qi(t) Qm 

qo{t) 



^^r^Mt) 



dq, 



i^ts {'^{'t^(l))^PsPs{v{t,q)) 



vr- 



■r^q'-qtit) 



go 



(t) JRe 



/^Ps {l{t,Q))T^PsPs{l{t,q)) 



IT- 



''rVqoit)-Q' 



dq 
dq, 



X V, 



Ps 



ifh<t<tj^ and- > , 

' ''max 

""PsPs.xV^'^'*) 

'^PsPs,y\^^'^^'^) 



' 12^ Ps 



90 



(t) 3f?e 



i7(t, q)K%^{-i{t, q))npsPs{l{t, q)) 



'Pt2Fts 



90 



(t) Ke 



/^Ps {l{t^(l))T^PsPs{l{t,q)) 



■dq. 



vrVVoP) 



dg, 



X y, 



X Vr 



if ti,_ < t and — > ^ or i/ to < ^ o'l-c^ ~ ^ 

' ^ max 



Ps 



and 



^ ''max 

upsPs{x, z,t) = e/se . 

VFe set /lere r = (x^ + (z + /i)^)-*^'^ and to = ''/^^s denotes the arrival time of th 
reflected PsPs volume wave at point (x,0,z), 



thi = (z + h) 



1 



1 \x\ 

+ 



Vpg ''max niiax 



(28) 



denotes the arrival time of the reflected PsPs head-wave at point {x, 0, z) and 



1 1 



Z + hVv+f ^max 



(29) 
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denotes the time after which there is no longer head wave at point {x, 0, z), (contrary to 
the 2D case, this time does not coincide with the arrival time of the volume wave). We 
also define the functions 7, v, qo and qi by 



xt z + h t"^ 1 

j:{teIR\t>to}xIR^C:= 7(t, Qy) = '^^ + ^-J ^ - ^JVT^^ 



V : {t e IR\thi <t < th2}>^IR^ C := v{t,qy) = -i 

qo : m ^ m := qo{t) 



z + h 



t"^ X 
+ ^t 



V^sily) "■' ' "■' 



\ 



t2 1 



r^ T/ 
y 1 



Ps 



and 



qi : IR ^ M := qi{t) 



\ 



-^\t-{z + h) 



V, 



+ 2 y2 



Ps 



y2 



ut, c is the solid displacement of the reflected PsS wave and satisfies: 



tip^5 3.(x,z,t) = — ^ / ^e 



^ ./O 



iv{t, q)K'^{v{t, q))npss{v{t, Q))-^(.t, q) 



dt 



dq, 



u'^^a-.{x,z,t) = — ^ / ^e 



-7r2 ; 
TT Jo 



dv 



{v\t,q) + q')npss{v{t,q))—{t, 



dq, 



■ifthi <t <to and 1 9m [7(^0,0)] I < 



^PsS.xi^'^'*) 



yxr 



+ 



'^Pss.zy^i^i'^) 



+ 



F^ 


rioit) 


•^Ps 

7r2 


/ ^e 
Jo 


T+ 


riiit) 


•^Ps 

7r2 


/ ^e 


Jqo{t) 


T+ 


riait) 


-' Ps 

7r2 


/ ^e 
Jo 


T+ 


m(t) 


'' Ps 

7r2 


/ ^e 



n{t,q)4{jit,q))npssil{t,q))^it,q) 
\v{t, q)K'^{v{t, q))npss{v{t, q))-K7{i, q) 



dt 



dq 

dq, 



dj 



{r{t,q) + qinpss{l{t,q))^{t,q) 
(v^t, q) + q^)npss{v{t, q))-Q^it, q) 



dq 
dq, 



if to < t < t/j2 o.'iT'd |9"i [7(^01 0)]| < 



V 



p+ j-qoit) 



-7r2 , 



i7(i, g)'ts' (7(i> q))'^pss{i{t, q))^{t, q) 



5t 



dq, 



p+ rqo{t) 
uj,^a-,{x,z,t) = — ^ / JRe 



7r2 ; 
^ JO 



97 



(7' (t, <?) + r )^Ps5(7(i, <?) ) ^ (t, (?) 



dg, 



INRIA 



Analytical solution 23 



"i-f ih2 < * ^'^'^ \Qm [7(^0,0)] I < — or if to < t and \Qm [7(^0) 0)] I ^ 



t' m a.x ►^ rr 

Pssy 



and Up a{x, z,t) = else. 



to denotes here the arrival time of the reflected PsS volume wave at point {x, 0, z) 
(its calculation is similar to the calculation of the arrival time of the transmitted wave, 
see the appendix of f^i), 



y V p max V 9 max ' max 

denotes the arrival time of the reflected PsS head wave at point (x, 0, z), 

,2 _L .2 _L ^^ /'c2 _L ci^ _L ^2 



h' + z' + hz{f + f\+x' 
*^. = n_\\ ^ (31) 

Cl C2 



denotes the time after which there is no longer head wave at point (2;,0, z), where 



Cl = , / ^ — 7-5 — and C2 



, F+/ ^^ax ' \ivf ^^ax' 

T/ie function qo : [to i +00] ^^ iR^ is the reciprocal function of to '■ IR^ >-^- [to, +00], 
where to{q) is the arrival time at point (x, 0, z) of the fictitious transmitted PsS volume 
wave, propagating at a velocity Vp^(g) from the source to the interface and at velocity 
Vg{q) from the interface to point {x,0,z) (we refer again to fS^ for details on its calcu- 
lation). 

The function qi : [ti ; to] ^^ ^^ is defined by 



Qiit) 



-^ I t- z 



\ '^ \ V S max V P s max / "^ ma 



The function 7 : {{t,q) £ IR^ x IR^ \ t > to{q)} ^^ C is implicitly defined as the only 
root of the function 

J-(7,g,t)=z(— i— + 72) +hl \ +^A +ijx-t 

whose real part is positive. 

The function v : Ei U E2 ^^ C is implicitly defined as the only root of the function 

( 1 \''' [ 1 \"' 

Tiv, q,t) = z\ T. Vv^\ + /i o Vv^\ + ivx — t 

such that Qm[dtv{t,q)] < 0, with 

El = {{t,q) £ IR+ X IR^\th, <t<to andO <q< qo{t)} 

and 

E2 = {(t,g) G ffi+ X iR+ I to < t < tft^ and qo{t) < q < qi{t)} . 
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Uprpr is the solid displacement of the transmitted PfPf wave (the Pf transm,itted 
wave generated by the Pf incident wave) and satisfies: 



^PfPfx\'^^ ^1 ) 



-PiiFpf P^W 



V/p/,2 



vr^ 



^e 



'x,z,t) 



Pn^P/ /■'^^(*) 



vr^ 



^e 



dv 

I ^^(i, q)'TpfPf{v{t, 9))-^(i, q) 



dv 



dq, 



Kpr{v{t, q))Tpfpf{v{t, q))jrr{t, q) 



dt 



ifthi <t ^to and 1 9m [7(^0,0)] | < 



yxr 



^PfPfxy'^^ ^1 ) 



KiFpf /•««(*) 



^PfPf^zy'^^ ^1 ) 



+ 



if to < t < th2 and |9m [7(^0, 0)]| < 



7r2 


Jo 


^ii^p/ 


/•gi(t) 


7r2 


Jqo(t) 


-PuF^f 


rio{t) 


7r2 


L 


'^iiFpf 


/•gi(t) 


7r2 


Jqoit) 


n^ll ^ 


1 



^e 



^e 



^e 



^e 



[-f{t,q)TpfPf{-/{t,q)) — {t,q) dq 

dv 
iv{t,q)Tpfpf{v{t,q)) — {t,q) dq, 

_ 97 

K-Pfhit, q))'FpfPf{i{t, 9))^(i, q) 

_ dv 

i^Pf(.y{t, q))'FpfPf{v{t, q))jrr{t, q) 



dt 



dq 
dq, 



•^rr 



''PfPf,x 



{x,z,t) 



KiFpf r*'W 



TT^ 



^e 



'PfPLz 



{x,z,t) 



KlFpf /■*'(*) 



vr^ Jo 



^e 



d^ 
n{t,q)TpfPf{-/{t,q)) — {t,q) 



d^ 



dq, 



i^Pfiiit, q))TpfPf{i{t, q))-Q^it, q) 



dq. 



ifth2 < * ^'^^ |9"i[7(*0,0)]| < — or if to < t and \Qm[j{to,0)]\ > — 

''max ''max 

anduprpAx,z,t) =0 else. 

to denotes here the arrival time of the transmitted PfPf volume wave at point (x, 0, z), 

(32) 



thi = h 



1 



1 



1 



F+Z K?.ax 



Z, 



1 \x\ 

+ 



Vp r ^max ^^max 



denotes the arrival time of the transmitted PfPf head wave at point {x,0,z), 



thn 



h _z_ 

Cl C2 



(33) 
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denotes the time after which there is no longer head wave at point {x,0,z), where 



1 



1 



ci 



V p r ' max 



and C2 



1 



1 



Vp r ' max 



The function qq : [tg ; +oo] h^ iR^ is i/ie reciprocal function of to : iR''' h^: [to,+oo], 
where tQ{q) is the arrival time at point (x, 0, z) of the fictitious transmitted PfPf volume 
wave, propagating at a velocity VpAq) in the top layer and at velocity VpAq) in the 
bottom layer. 

The function qi : \ti ; to] ^^ ^^ ^-^ defined by 

«i(«) 



\ 



^ t + z. 



V p r ' max 



V+; ^^ax 



y2 

' ma 



The function 7 : {(t,g) G iR+ x iR+ 1 1 > to(g)} 
root 0/ t/ie function 



is implicitly defined as the only 



1/2 



^h,<l,t) = -z 



Vp/{q) 



+ 7^ 



+ h 



1 2 



1/2 



V^/(g) 



+ ijx — t 



whose real part is positive. 

The function v : EiD E2 ^^ C is implicitly defined as the only root of the function 



J='{v,q,t) 



1 



1/2 



Vp/{q) 



+ v^ 



+ h 



1 



1/2 



.V+/(g) 



+ v' 



+ ivx — t 



such that Qm [(9it; (t, g)] < 0, with 

El = {{t,q) € m^ X lR+\th^<t<to andO<q< qo{t)} 

and 

E2 = {it,q) £ IR'^ X IR'^\to <t <th^ and qo{t) < q < qi{t)} . 

Uprpg is the solid displacement of the transmitted PfPs wave and satisfies: 



^PfPs,a;V^'^'*) 



Vi^F+j nAt) 



vr^ 



^e 



u 



PfPs,z 



ix,z,t) 



'PuFpf /"^^(*) 



vr^ 



3f?e 



dv 
iv{t, q)TpfPs{v{t, Q))-g^{i, q) 



dv 



dq, 



'^Ps(^(*' QWpfPs{v{t, Q))-Q^{t, q) 



dq, 
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ifthi <t<to and 1 9m [7(^0,0)] | < 



V ' 

y max 



U 



PfPs,x 



[x,z,t) 



Vi2Fpf /■*'(*) 



vr^ 



^e 



v^^Fpf r"^^'^ 



vr^ 



^e 



90 (t) 



U 



PfPs,z 



'x,z,t) 



VuFpf /•*'(*) 



vr^ 



^e 



+ 



^r2i^P/ /"^^w 



^e 



iftQ<t< th2 and l^m [7(^0,0)]! < 



1 



57 
['y{t,q)TpfPs{'y{t,q)) — {t,q) dq 

dv 
iv{t,q)TpfPs{v{t,q)) — {t,q) dq, 



'^Ps(7(i, qWpfPs{i{t, 9))-^(i, q) 



_ dv 

i^Psi^it, q))TpfPs{v{t, Q))-Q^{i, q) 



dq 
dq, 



V 

•'max 



'^PfPs,x\^^^^^) 



V^2Fpf /■^«(*) 



vr^ 



Ke 



'^PfPs,z\^-''^i'^) 



K^Fpf /"^»(*) 



vr^ Jo 



^e 



i7(t, q)TpfPs{-f{t, Q))-g^{i, q) 



dj 



dq, 



i^Psiiit^ (iWpfPsiiit, 9))-^(i, q) 



dq. 



ifth2 < * ^'^^ |9"i[7(*o,0)]| < — — or ifto<t and \Qm[y{to,0)]\ > — — 

''max ''max 

and Uprp^{x,z,t) = else. 

to denotes here the arrival time of the transmitted PfPs volume wave at point (x,0, z), 



thi = h 



1 1 



y+/ ^rLx 



1 \x\ 

+ 



Vpg Viax ''max 



(34) 



denotes the arrival time of the transmitted PfPs head wave at point (x,0, z), 



th2 



K' + z'-hz[f^+f^]+x' 

h _z_ 

Cl C2 



(35) 



denotes the time after which there is no longer head wave at point (2;, 0,2;), where 



Cl 



y+/ ^^ax 



and C2 



1 1 



V pg ^ max 



The function qo : [to; +00] 1-^ IR^ is the reciprocal function of to : IR^ 1-^: [to, +00], 
where to{q) is the arrival time at point (x, 0, z) of the fictitious transmitted PfPs volume 
wave, propagating at a velocity Vpr{q) in the top layer and at velocity Vp^(g) in the 
bottom layer. 
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The function qi : [ti ; to] ^^ ^^ ^-5 defined by 

m{t) 



\ 



1 1 1 

t + Z 



1 1 



V n „ "^ max \ V m III 



max / "^ max 



T/ie function 7 : {(t,Q') G iR''' x IR^ \ t > ^0(9)} ^^ C is implicitly defined as the only 
root of the function 



J^ii,q,t) = -z 



1 



1/2 



1 



1/2 



, +7 + ^ 



2 , , + 7 I + hx - 1 



pf 



whose real part is positive. 

The function v : Ei U E2 ^^ C is implicitly defined as the only root of the function 



J^{v,q,t) = -z 



1 



Vp/(g) 



\ 1/2 



1 



.V+/(g) 



2 V'' 

+ V \ + ivx — t 



such that 9m[5if(t, g)] < 0, with 

El = {{t,q) e m^ X lR^\th^<t<to andO < q < qo{t)} 

and 

E2 = {it,q) £ R^ X M^ \to < t < th^ and qo{t) < q < qiit)} . 

Up fa is the solid displacement of the transmitted PfS wave and satisfies: 



'■PfS 



^PfS,x(^'^'*) 



F+^ r^iW 



^PfS,z\^^ ^' / 



^ Jo 






^e 



_ dv 

yv{t, q)Kg {v{t, q))Tpfs{v{t, Q))-g^it, q) 



dq, 



^e 



{v\t,q) + q^)Tpfsiv{t,q)) — {t,, 



dq, 



ifthi <t <to and 1 9m [7(^0,0)] | < 



Vrr- 



^P/S,x(^'^'*) 



F+^ r*)(t) 



-7r2 , 
^ ./O 



^e 



Ftr rii{t) 



7r2 



^e 



go{t) 



_ ^7 

n{t, q)Ks hit, q))'Fpf shit, q))—{t,q) 

_ dv 

ivit, q)Kg ivit, q))Tpfsivit, q)) — it, q) 



dt 



dq 
dq. 



'^PfS,zi-^i Z,t) 



T^j r^oct) 



-7r2 , 



^e 



dj 



+ 



T^f r^iW 



vr^ 



3fie 



-?(){*) 



iYit,q) + q')Tpfsilit,q))^it,q) 
iv\t, q) + q')rpfsivit, q))^it, q) 



dq 

dq. 
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if tQ < t < t/j2 and \^m[y{to,0)]\ < 



yxr 



U 



PfS,x 



{x,z,t) 



1_P/ 

-7r2 , 
^ ./O 



90 (t) 



Re 



n{t, q)Kg (7(t, q))Tpfs{j{t, q))-K7{t, q) 



dt 



dq, 



^PfS.zy^^ ^' ^) 



T\ 



Pf 



^ ./O 



90 (t) 



Re 



dj 



{Y{t,q) + q')Tpfs{l{t,q))^{t,q) 



dq, 



if ih2 < * ^'^^ I^^T- [7(^0) 0)]| < — or if tQ <t and |9m [7(^0, 0)]| > 



Vrn 



v„. 



and Uprg{x, z, t) = else. 

to denotes here the arrival time of the transmitted PfS volume wave at point {x,0,z), 

(36) 



thi 



h. 



1 



1 



1 



V+; K^ax 



1 \x\ 

+ 



Vs' K?.ax ' K.ax 



denotes the arrival time of the transmitted PfS head wave at point (x,0,z), 

th2 



h'+z'-hz{^ + ^]+x' 



h _z_ 



(37) 



denotes the time after which there is no longer head wave at point (x,0, z), where 



1 



1 



Cl 



v+; ^4 



and C2 



1 



1 



V. 



_2 



1/2 



The function go '■ [to] +00] h^ IR^ is the reciprocal function of Iq : IR^ h^: [to, +00], 
where to^q) is the arrival time at point (x, 0, z) of the fictitious transmitted PfS volume 
wave, propagating at a velocity VpJq) in the top layer and at velocity Vg{q) in the 
bottom layer. 

The function qi : [ti ; to] ^^ ^^ is defined by 
9l(«) 



\ 



"* + .," 



1 



K 



_2 



1/2 



T/+ 2 T/2 
V p r ' max 



V2' 



The function 7 : {{t,q) G SV' x iR^ 1 1 > to(g)} ^^ C is implicitly defined as the only 
root of the function 



^ii,q,t) 



1 



1/2 



ys (q) 



+7^ 



+/i 



1 



1/2 



v^/(<z) 



+ 7^ 



+ i7x — t 



whose real part is positive. 

The function v : Ei U E2 ^^ C is implicitly defined as the only root of the function 



F{v,q,t) = -z 



yfi<i) 



1/2 



+ v' 



+ h 



1 



1/2 



+ v' 



+ ivx — t 
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such that Qm [dtv{t, q)] < 0, with 

El = {{t,q) e m^ X lR+\th, <t<to andO <q< qo{t)} 

and 

E2 = {{t,q) e M'^ X m^ \to < t < th^ and qo{t) < q < qi{t)} . 

Up^pr is the solid displacement of the transmitted PsPf wave and satisfies: 






vr^ Jo 



dv 
i v{t, q)TpsPfiv{t, Q))-g^{t, q) 



dq, 



'^PsPf,z\^^^^^) 



-PiiFps r^'^ 



vr^ 



^e 



_ dv 

i^Pfi^^ii, q))'FpsPf{v{t, 9))iTr(i, 



dt 



dq, 



ifthi <t <to and \Qm[y{to,0)]\ < 



Kr 



U 



PsPf,x 



'x,z,t) 



-PuFps r^'^ 



^e 



TT- JO 



vr^ 



^PsPf.^i^'^'*) 



'go(t) 



+ 



T^ Jo 



vr^ 



^e 



^e 



^e 



i7(t, q)TpsPf{-/{t, Q))-g^{t, q) 



dv 
iv{t,q)TpsPf{v{t,q))—-{t 



dt 



dq 
dq, 



dj 



Up f (lit, qWpsPfiiit, Q))-K7it, q) 



if to < t < t/i2 and \Qm [7(^9, 0)]| < 



90 (t) 
1 



dv 



KpAv{t, q))TpsPf{v{t, q))-^{t, q) 



dt 



dq 
dq. 



V ' 



UpsPf.ix,z,t) = '-^^ I ^e 



vr^ 



^PsPf^zy^i'^-''^) 



KiFps r^'^ 



vr^ 



^e 



d^ 
\-f{t,q)TpsPf{-i{t,q)) — {t,q) dq, 

_ d^ 

i^Pf{i{t, q))TpsPf{-t{t, q))^{t^ q) 



dq. 



if th2 < t and \Qm [7(^0, 0)]| < — or if t^ <t and |9m [7(to,0)]| > 



Vrn 



v„ 



and Up^pAx,z,t) = else. 

to denotes here the arrival time of the transmitted PsPf volume wave at point (x, 0, z) 



thi — h. 



V, 



Ps 



2 y2 



1 \x\ 

+ 



Vp r max "^ max 



(38) 
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denotes the arrival time of the transmitted PsPf head wave at point (x, 0, z), 



thn 



h _z_ 

ci C2 



(39) 



denotes the time after which there is no longer head wave at point {x,0,z), where 



ci 



Vd„ ''^ma 



and C2 



V, 



_ 2 



1/2 



The function qq : [to; +oo] h^ IR^ is the reciprocal function of to : IR^ ^^: [to,+oo], 
where to^q) is the arrival time at point (x, 0, z) of the fictitious transmitted PsPf volume 
wave, propagating at a velocity Vp^{q) in the top layer and at velocity VpAq) in the 
bottom layer. 

The function qi : \ti ; to] ^^ ^^ ^-^ defined by 



9l(() 



\ 



-^ t + z. 



v^; ^.Lx 



V p^ ^ max 



y2 • 



The function 7 : {{t,q) G IR^ x IR^ \ t > tQ{q)} h^ C is implicitly defined as the only 
root of the function 



^{i,q,t) 



1 



Vp/{q) 



+ 7 



\ 1/2 
2^ +h 



1 



ykiQ) 



2 V'' 

+ 7 + ^7X — t 



whose real part is positive. 

The function v : Ei L) E2 ^^ C is implicitly defined as the only root of the function 



J='{v,q,t) = -z 



1 



1/2 



Vp/(g) 



+ v' 



+ h 



1 



1/2 



yts\<i) 



+ v' 



+ ivx — t 



such that Qm[dtv{t,q)] < 0, with 

El = {{t,q) e m^ X lR+\th, <t<to andO < q < qo{t)} 

and 

E2 = {it,q) £ IR^ X IR^\to <t <th^ and qo{t) < q < Qiit)} ■ 

• Upgpg is the solid displacement of the transmitted PsPs wave and satisfies: 



""PsPs.a:!^'^'*) 



■p- /?+ fjiW 



■K" 







dv 
\v{t, q)TpsPs{v{t, (l))-Q^{t, 



dq, 






vr^ 



^e 



_ dv 

i^Ps(.'^ii^'i))'^PsPs{v{t,q)) — {t, 



dq, 
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ifthi <t<to and 1 9m [7(^0,0)] | < 



^PsPs,a;(^'^'*) 



'^PsPs.zy^i'^i'^) 



+ 



if tQ < t < th2 and \Qm[y{tQ,0)]\ < 



h^Jil ^ 


''max 


7r2 


f-qoit) 
JO 


^12 -^Ps 

7r2 


Jqo{t) 


7r2 


rqo(t) 

^0 


n-2^i 

7r2 


Jqo{t) 


n^ll ^ 


1 



^l{t,q)TpsPs{j{t,q))-^{t,q) dq 

dv 
iv{t,q)TpsPs{v{t,q)) — {t,q) dq, 

_ d^ 

i^Pshi^^ qWpsPs{i{t, 9))-^(i, q) 



_ dv 

i^Psi'^i.i^ q))TpsPs{v{t, Q))-Q^{t, q) 



dq 
dq, 



V 



V^ /?+ rio{t) 



vr^ 



-p^ p+ f-qoit) 



vr^ Jo 



i7(t, q)TpsPs{-i{t, Q))-Q^(.i^ ?) 



97 



dg, 



/^PsItI^, q))TpsPs{l{t, Q))-g^it^ Q) 



dq. 



ifthi < * '^'^d 1 9m [7(^0,0)] I < — or if to < t and 1 9m [7(^0,0)] | > — 

''inax ''max 

and Up^p^{x, z,t) = else. 

to denotes here the arrival time of the transmitted PsPs volume wave at point (x, 0, z), 

(40) 



thi = h. 



1 



Vp\' "^ 



1 1 X 

+ 



max \/ Vr,„ "^nax ''max 



Ps 



denotes the arrival time of the transmitted PsPs head wave at point {x,0,z), 



thn 



h _z_ 

Cl C2 



(41) 



denotes the time after which there is no longer head wave at point {x,0,z), where 



1 1 



Cl 






and C2 



1 



Ps 



1/d„ "^max 



Ps 



The function qo : [to] +00] h^ IR^ is the reciprocal function of to : IR^ ^^: [toj+oo], 
where to{q) is the arrival time at point (x, 0, z) of the fictitious transmitted PsPs volume 
wave, propagating at a velocity Vpg{q) in the top layer and at velocity Vp^(g) in the 
bottom layer. 
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The function qi : [ti ; to] ^^ ^^ ^-5 defined by 

«i(«) 



\ 



1 / 1 






' -/..M 



'Ps 






1/2 • 

* vn a-v- 



T/ie function 7 : {(t, Q') € IR^ x iR^ 1 1 > tQ{q)} h^ C is implicitly defined as the only 
root of the function 



J^ii,q,t) = -z 



1 



,- 2 



VpsiQ) 



\ 1/2 

+ 7M +h 



1 



v^'(g) 



2 V'' 
+ 7 + i^x — t 



whose real part is positive. 

The function v : Ei U E2 ^^ C is implicitly defined as the only root of the function 



J='{v,q,t) = -z 



1 



.Vp/(g) 



\ 1/2 
+ v^] +h 



1 



yk'iQ) 



2 V'' 

+ V \ + ivx — t 



such that Qm [dtv{t, q)] < 0, with 

El = {{t,q) e IR+ X m+\th, <t<to andO <q< qo{t)} 

and 

E2 = {it,q) G -K+ X iR+ I to < t < t/ii and qo{t) < q < gi(t)} . 

• Up^g is the solid displacement of the transmitted PsS wave and satisfies: 



Fp /■9i(*) 



VT ./o 



_ dv 

iv{t, q)Kg {v{t, q))Tpss{v{t, 9))-^(*> O) 



dq, 



Up^a-.{x,z,t) = ^TT- I ^e 



-7r2 , 
^ ./O 



(vHt, q) + q^)Tpssivit, q))^it, q) 



dq, 



ifthi <t <to and |3=m [7(to,0)]| < 



^Ps5,x(^'^'*) 



Kr 



^Ps5.^(^'^'*) 



+ 



7r2 


/ ^e 
Jo 


7r2 


riiit) r 
/ ^e 


7r2 


Jo 


■^Ps 

7r2 


/•9l(*) 
Jqoit) 



n{i, q)i^s (7(*> qWpssiiit, ?))"q^(*' «) 



_ 9f 

if (t, g)K_5 (i;(t, q))Tpss{v{t, q))^{t^ q) 



dq 
dq, 



{l^{t,q) + q^)Tpss{l{t,q))^{t,q) 
{v\t, q) + q^)Tpss{v{t, q))^it, q) 



dq 
dq. 
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if tQ < t < t/j2 and \Qm[y{to,0)]\ < 



yxr 



^Ps5.x(^'^'*) 






-7r2 , 



n{i, q)i^s (7(*> QWpsshit, 9))^{t, q) 



dt 



dq, 






^ Jo 



d^ 



{rit,q) + q')rpss{l{t,q))-^{t,q) 



dq, 



ifth2 <* o.^d |9m[7(to,0)]| < — oriftQ <t and \Qm[y{to,0)]\ > — 

''max ''max 

and Up_,g{x, z,t) = else. 

to denotes here the arrival time of the transmitted PsS volume wave at point {x,0,z), 

(42) 



thi 



h. 



1 



1 



1 



^Ps ^°^« 



1 \x\ 

+ 



Vs^ ^^ax ' ^max 



denotes the arrival time of the transmitted PsS head wave at point (x, 0, z) 



tho 



h^ + z^-hz(^ + ^] +x2 

h _z_ 

Cl C2 



(43) 



denotes the time after which there is no longer head wave at point (x, 0, z), where 



1 



1 



Cl 



V, 



+ 2 y2 



and C2 



1 



1 



Ps 



V, 



-2 y2 



The function qo : [to; +oo] >-^ IR^ is the reciprocal function of Iq : IR^ h^: [to,+oo], 
where to{q) is the arrival time at point {x, 0, z) of the fictitious transmitted PsS volume 
wave, propagating at a velocity Vp^(g) in the top layer and at velocity Vg (q) in the 
bottom layer. 

The function qi : \ti ; to] ^^ J^^ is defined by 



ii(t) 



\ 



-^ { t + z, 



1 



1 



1 



T/-2 y2 

V Q ' ma 



The function 7 : {{t,q) £ IR+ x IR+ \t > io{q)} 
root of the function 



T/+2 T/2 I y2 • 

V D „ max / "^ max 

is implicitly defined as the only 



^{i,q,t) = -z 
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whose real part is positive. 

The function v : EiD E2 ^^ C is implicitly defined as the only root of the function 



F{v,q,t) = -z 



1 



1/2 



Vc (9) 
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such that Qm [dtv{t, q)] < 0, with 

El = {{t,q) e m^ X lR+\th, <t<to andO <q< qo{t)} 

and 

E2 = {{t,q) e M'^ X m^ \to < t < th^ and qo{t) < q < qi{t)} . 

Remark 2.1. For the practical computations of the velocities, we won't have to explicitly 
compute the derivatives of the displacement u, which would be rather tedious, since 

du , 

dt •' •' 

Therefore, we'll only have to compute the derivative of the source function f. 

3 Numerical illustration 

To illustrate our results, we have computed the green function and the analytical solution to 
the following problem: we consider an two-layered poroelastic medium whose characteristic 
coefficients are 

• the solid density: p^ = 2200 kg/m^ and p^ = 2650 kg/m^; 

• the fluid density: pt = 950 'kg/va? and p7 = 750 kg/m^ ; 

• the porosity: 0^ = 0.4 and (j)~ = 0.2 ; 

• the tortuosity: a^ = 2 and a^ = 2; 

• the solid bulk modulus: Kf = 6.9 GPa and Kj = 37 GPa; 

• the fluid bulk modulus: KJ = 2 GPa and KJ = 1.7 GPa; 

• the frame bulk modulus: K^ = 6.7 GPa and K^ = 2.2 GPa; 

• the frame shear modulus p+ = 3 GPa and p," = 4.4 GPa; 
so that the celerity of the waves in the poroelastic medium are: 

• for the fast P wave, Vp^ = 2692 m/s and Vpp = 2535 m/s; 

• for the slow P wave, Vp^ = 1186 m/s and Vp^ = 744 m/s; 

• for the ifj wave, Vg = 1409 m/s and Vg = 1415 m/s. 

The source is located in the top layer, at 500 m from the interface. We used two types of 
sources in space: the first one is a bulk source such that /« = /«, = — lO^'^ and fp = 0; the 
second one is a pressure source such that fu = fw = and /p = 1. In each case we used a 
fourth derivative of a Gaussian of dominant frequency /q = 15 Hz: 



7r2 



m = 272 



2/ i\2 4/ i\4' 

vr / 1 \ . TT / 1 ^ 



3 + 12-2 t-- +4-4 t-- 



'l^y ft / 
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for the source in time. We compute the solution at two receivers, the first one is in the upper 
layer, at 533 m from the interface; the second one is in the bottom layer, at 533 m from 
the interface; both are located on a vertical line at 400 m from the source (see Fig. [3]). We 
represent the z component of the green function associated to the solid displacement from 
t = to i = 1.4 s in Fig. m for the bulk source and in Fig. E] for the pressure source. In Figs. 
andm we plot the solid displacement. The left pictures represents the solution at receiver 1 
while the right pictures represents the solution at receiver 2. As all the types of waves are 
computed independently, it is easy to distinguish all of them, as it is indicated in the figures, 
solution. 



Receiver 1 



n+ 



Source 
500 m 



400 m 



533 m 



Q- 



533 m 
Receiver 2 



Figure 3: Configuration of the experiment 
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Figure 4: The z component of the green function associated to the displacement at receiver 
1 (left picture) and 2 (right picture), in the case of a bulk source 
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Figure 5: The z component of the displacement at receiver 1 (left picture) and 2 (right picture) 
in the case of a bulk source. 
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Figure 6: The z component of the green function associated to the displacement at receiver 
1 (left picture) and 2 (right picture) in the case of a pressure source. 
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Figure 7: The z component of the displacement at receiver 1 (left picture) and 2 (right picture) 
in the case of a pressure source. 
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